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KEY WORDS: Alaska; Fjords; Inlets (waterways); Ocean currents; Ocean 
engineering; Oceanography; Water transfer 


ABSTRACT: part comprehensive program marine environmental studies 
Port Valdez, Alaska, currents were measured near Valdez Narrows almost 
continuously for more than two years. The purpose these measurements was 
investigate the exchange water between Port Valdez and Prince William Sound. 
Intensive profiling temperature and salinity Port Valdez and adjacent waters has 
enabled observation water exchange processes that are effective seasonal basis. 
However, the current meter data have shown that these processes are occasionally 
overshadowed short-term events that are comparable effectiveness promoting 
the exchange water between Port Valdez and Prince William Sound. The latter 
events have typical durations two five days and their timing apparently related 
the passage weather systems through the Gulf Alaska. The random nature 
these events, conjunction with the generally predictable seasonal exchange 
processes, suggests that the notion unique residence time for Port Valdez and 
similar fjords inappropriate. 


REFERENCE: Colonell, Joseph (Senior Project Woodward-Clyde 
Consultants, Anchorage, Alaska 99503), Water Renewal Subarctic Fjord,” 
Journal the Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 107, No. 
Proc. Paper 16632, November, 1981, pp. 223-231 
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KEY WORDS: (waterways); Models; 
Sedimentation; Shoals; Tidal traffic flow; Tides 


ABSTRACT: The morphology ebb tidal shoals has received considerable attention 
the last few years, but generally applicable techniques have evolved predict 
shoal development. This engineering study first step that direction. Shoals were 
generated laboratory basin and the basin controlling parameters were identified. 
Results indicate that the inlet behaves two dimensional turbulent jet and that 
inlet width, discharge velocity, and critical sediment velocity govern the equilibrium 
shoal dimensions. This preliminary study also suggests that the development shoal 
length and width are governed different transport mechanisms. Resolution such 
possibilities requires more detailed studies. 


REFERENCE: Sill, Ben (Assoc. Prof., Dept. Civ. Engrg., Clemson Univ., 
Clemson, C.), Fisher, John S., and Whiteside, Stuart D., “Laboratory Investigation 
Ebb Tidal Shoals,” Journal the Waterway, Port, Coastal and Ocean 
ASCE, Vol. 107, No. WW4, Proc. Paper 16674, November, 1981, pp. 233-242 
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WATER RENEWAL SUBARCTIC FJORD 


Port Valdez, Alaska, has gained prominence the location the marine 
terminal the Trans-Alaska Pipeline System. part acomprehensive program 
marine environmental studies this subarctic fjord, currents were measured 
Valdez Narrows for more than two years examine the exchange water 
between Port Valdez and Prince William Sound. These measurements have shown 
that flushing Port Valdez probably both frequent and intensive but the 
notion unique residence time for the fjord inappropriate. 

Physical Oceanography Port Valdez.—Port Valdez and Valdez Arm form 
(46 km) northeasterly extension Prince William Sound. Communication 
with the north Pacific Ocean via Hinchinbrook Entrance, which opens 
the Gulf Alaska about (110 km) from the town Valdez (Fig. 1). 
The fjord Port Valdez extends about (23 km) eastward from Valdez 
Narrows, narrow, double-silled entrance that regarded the demarcation 
between Port Valdez and Valdez Arm. The town Valdez lies near the northeast 
corner the fjord while the Trans-Alaska Pipeline System marine terminal 
the south shore Jackson Point (Fig. 2). 

Shaped somewhat like bathtub, Port Valdez about km) wide 
(18 km) long with steep sides the north and south and nearly 
horizontal bottom depth about 790 (240 over three-quarters 
its length. the easternmost quarter the fjord, the bottom rises rather uniformly 
the eastern shore the former townsite Valdez. The maximum depth 
Port Valdez 810 (247 m), while the overall mean depth the fjord 
about 580 (180 m). 

Port Valdez has oceanographic regime that strongly stratified (vertically) 
both temperature and salinity during the summer months but virtually 
unstratified during the winter months due thermohaline convection (2,6). 
This fjord would classified estuary since receives more 
fresh water runoff and precipitation than lost through evaporation (1,2). 
This should imply that the classical estuarine circulation, with seaward movement 
brackish surface layer and landward movement deeper waters, would 


Project Engr., Woodward-Clyde Consultants, Anchorage, Alaska. 

Note.—Discussion open until April 1982. extend the closing date one month, 
written request must filed with the Manager Technicai and Professional Publications, 
ASCE. Manuscript was for review for possible publication August 1980. 
This paper part the Journal the Waterway, Port, Coastal and Ocean Division, 
Proceedings the American Society Civil Engineers, Vol. 107, No. WW4, 
November, 1981. ISSN 0148-9895 /81 /0004-0223 $01.00. 
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prevail Port Valdez least during periods maximum runoff. Hydrographic 
data suggest that such flow does occur during the summer months but that 
confined the top (15 the water column (2,6). 

The tides Port Valdez are the mixed, semidiurnal type with maximum 
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FIG. 2.—Map Port 


range 17.5 (5.3 and mean height 10.0 (3.0 m). The tidal prism 
about 1.6% the total volume Port Valdez (2). Assuming homogeneous, 
well-mixed water body, i.e., winter conditions, can calculated that the 
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half-life volume conservative substance about days Port Valdez. 
This ignores any benefit dilution fresh water input, which would negligible 
under winter conditions. 

Just south Valdez Narrows there are two sills, with maximum depth 
525 (160 that limit direct exchange Port Valdez water below that depth 
with the deep waters Prince William Sound. should also noted that 
direct exchange water between the Sound and the Gulf Alaska similarly 
inhibited sill with maximum depth about 574 (175 south 


a 


(m) 


DEPTH 


FIG. 3.—Temporal Variation Water Density Hydrographic Stations: (a) 
Valdez Arm (VAA, upper); (b) Port (PVA 139, lower) (Specific Gravity 
10° 


Hinchinbrook Entrance. Temporal variations water density for one station 
within Port Valdez (PVA 139) and one Valdez Arm (VAA 110) are presented 
Fig. These figures are essentially maps the variation water density 
with time. Each vertical dashed line corresponds the density profile observed 
the indicated time. The contours have been drawn simply connecting 
points equal density without attempting smooth the data for more realistic 
estimate the variations between observations. Also, the difference vertical 
scales should noted. 
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Evidence the intrusion dense water from Prince William Sound into 
Valdez Arm and Port Valdez apparent Fig. Appearing pulses 
increased density mid-to-late summer, these annual intrusions serve replenish 
bottom waters several fjords southcentral Alaska (5). The driving force 
for these intrusions has been attributed the combined effects prevailing 
winds over the Gulf Alaska and the voluminous precipitation-runoff the 
southcentral Alaska coastal region (8,9). During the summer and early autumn 
months, deep water exchange enhanced the development classical 
estuarine flow conditions. The individual effectiveness these processes 
difficult determine because their simultaneous occurrence but prior study 


TABLE 1.—Valdez Narrows Current Meter 


Current Speeds” 


Stan- 

dard 
devia- Maxi- 

Prin- Mean, tion, 

cipal 
Water Meter Useful axis centi- centi- centi- 


(3) (4) (5) (6) (7) (8) (9) 


1977 
April 19, 
1978 

July 16, 
1978 
November 
11, 1978 


moorings were located 0.6 (0.9 km) Middle Rock (61.07° 146.66° W). 
Aanderaa RCM-4 current meters were used with data sampling interval set min. 

from unfiltered data. VN-4 results are for data records truncated 19-day 
period during which all four current meters were operating (December, 8-26, 1977) 

“Short records due battery failures. 


Port Valdez suggested residence times approx days (6). 

Valdez Narrows Currents.—During the period from October, 1976-November, 
1978, four long-term current meter moorings were placed near Valdez 
Narrows. The purpose these moorings was measure the exchange water 
between Port Valdez and Prince William Sound. Although data from all four 
moorings have been examined, the most interesting and useful information has 
been derived from the two moorings that were deployed the sill Valdez 
Narrows (VN and #5) the location shown Fig. Some operational 
details for these moorings are provided Table 

Current meter data were analyzed, first, terms true north-south and 
east-west components the velocity vector. For data that are obtained 


Moor-| Dates 
ing 
deploy- 
ber ment 
106 132 052 
156 34° 035 
110 118 061 
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well-defined channel, such Valdez Narrows, generally more useful 
have the velocity vector resolved into its along-channel and cross-channel 
components. Rather than assume channel orientation from its appearance 
nautical chart, the statistics the velocity data are used obtain the preferred 
direction flow, which termed the the flow (4). The 
latter defined mathematically that axis orientation for which the covariance 
the velocity components minimized. shown Table the principal 
flow directions for the current meters single mooring tended quite 


SPEED (cm/sec) 


CURRENT 


° 


i977 


FIG. 4.—Along-Channel Velocity Component Narrows (VN#4, 10-28 Dec., 
1977). Positive Speed Northeasterly Flow (Azimuth 46°); Negative Speed 
Southwesterly Flow (Azimuth 226°). Signal Shown Fine Line; 
Dark Line Net Speed after Removal Tidal Constituents Digital Filtering 


similar, i.e. within few degrees, that average angle for coordinate axis 
rotation was used obtain the along-channel and cross-channel velocity compo- 
nents (termed and respectively) for all meters. Some statistics observed 
current speeds are also listed Table All subsequent computations were 
performed the velocity components obtained for the principal axis 
orientation. 


obtain information the net exchange water through Valdez Narrows, 
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was necessary remove the tidal signal from the raw current meter data. 
Butterworth low-pass filter was used for this purpose (7) with example 
the results (Fig. 4). the latter figure both filtered and unfiltered data 
are displayed for the along-channel velocity component each four depths 
obtained during December, 1977 (VN #4), showing how the mean flow 
extracted from the total flow signal. 

Because battery failures the current meters used the mooring, 


(cm/sec) 


SPEED 


WwW 
a 
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FIG. 5.—Along-Channel Velocity Component Valdez Narrows (VN#5, 
November 1978). Positive Speed Northeasterly Flow (Azimuth 50°); Negative 
Speed Southwesterly Flow (Azimuth 230°) 


the period during which all four meters were operating was limited December 
10-28, 1977. indicated Table useful data were obtained from other 
meters this mooring for longer period; however, sufficient for this 
discussion focus these shorter records. Referring Fig. two features 
are noteworthy. First, classical estuarine flow apparent during the period 
December 10-17, despite the fact that fresh water input virtually zero 
that time year. Second, for the remainder the record, the mean flow 


= 
| 28 m 
20 ai 
40 
20 66 m 
-2 4 
V 
-40: 
2 
~40: 
' 
20 152 m 
-40 
| 
At 


ww4 DEEP WATER RENEWAL 229 


profile reversed with surface layer flow into Port Valdez and outward 
the deeper layers. Weather data, consisting winds Middle Rock and daily 
mean barometric pressure Valdez, were compared with the mean flow variations 
that are indicated Fig. was concluded that the latter are not directly 
related the local winds measured Middle Rock, but are fairly responsive 
barometric pressure This responsiveness was especially prominent 
for the period December 12-16, when two-lobed variation current speed 
all four meters coincided with similar variation barometric pressure 
(2). 

Fig. the along-channel velocity components are plotted for all four current 
meters mooring for the period July 12, 1978. While 


ALONG-CHANNEL CURRENT SPEEDS (cm/sec) 


Southwesterly Northeasterly 
3 (out of Port Valdez) (into Port Valdez) 
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FIG. 6.—Cross-Section Valdez Narrows, Showing Area Allocation Current Meters 
for Flow Computations (Left); and Velocity Profiles (Net Speeds) 
Observed Indicated Dates (Right) 


occurrences classical estuarine flow were definitely evidence, would 
expected that time year, the prominent features these records are 
the opposite flow conditions; that is, there were numerous occurrences flow 
into Port Valdez the surface layer with simultaneous outward flow the 
deeper layers. While the September 12-14 cannot escape 
notice, should noted that less spectacular but similar events occurred 
several times during the observation period. Attempts correlate such events 
with available weather data were less satisfying than for the previous example 
and had regarded inconclusive. Nevertheless, the hypothesis causal 
relationship between these events and the local regional weather conditions 
reasonable and warrants additional study. For purposes this discussion, 
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however, sufficient recognize the occurrence these events and their 
probable relationship with weather patterns that are typical southcentral Alaska. 
remains estimate the role these events promoting the exchange 
water between Port Valdez and Prince William Sound. 

Fig. includes three examples along-channel velocity profiles indicative 
surface inflow-deep outflow events observed during December, 1977, and 
September and October 1978. the basis these profiles, was assumed 
that the depth zero net motion was about should emphasized 
that the indicated shapes these velocity profiles are only estimates based 
experience with variety flow phenomena; however, the profiles are 
probably satisfactory for the outflowing layer where the flow reasonably 
well defined observations three depths. was further assumed that depths 
and 140 were reasonable limits for the channel cross-sectional areas 
associated with the three deepest current meters mooring #5, 
least for these flow events. 

Using these allocations channel cross-sectional areas (Fig. and some 
fairly conservative estimates mean flow speeds and flow durations, the net 
outflow water from Port Valdez through Valdez Narrows was calculated 


TABLE 2.—Computed Outflows from Port 


Total outflow 
(10° 
(2) 


Outflows percentage 
Port Valdez volume 


December 20-28, 1977 
September 13-15, 1978 
October 6-8, 1978 


for several the observed the current meter records. Computations 
for the period December, 1977 (from mooring #4) were essentially the 
same for September and October, 1978 except that the depth zero net 
motion was assumed the second current meter (Fig. 4). Examples 
the computed outflows are listed Table With magnitudes exceeding 
15-20% the volume the fjord, apparent that only few such events 
would required flush Port Valdez rather thoroughly. 


Surges surface water from Prince William Sound into Port Valdez, ac- 
companied large outflows from depths greater than 50-65 appear 
related the passage weather systems through southcentral Alaska. These 
events appear more effective promoting deep water exchange than 
either tidal action annual bottom water intrusions. Because their apparent 
connection with atmospheric conditions, certain degree randomness 
associated with these events that characteristic flushing rate difficult 
identify; however, the latter probably matter more than several 
days, view the large outflow volumes that were indicated the current 
meter data. 
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The random nature the weather-related events, conjunction with the 
generally predictable seasonal and tidal exchange processes, implies that the 
notion unique residence time for Port Valdez inappropriate. However, 
the observational evidence suggests that conservative estimate composite 
residence time probably more than few weeks. The latter depends, 
course, assumption reasonably uniform mixing throughout Port 
Valdez. 

These conclusions should apply equally well other fjords southcentral 
Alaska, and perhaps some southeastern Alaska and northern British 


Columbia, since most are subject the same meteorologic and oceanographic 
forcing functions. 
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LABORATORY INVESTIGATION EBB TIDAL SHOALS 


Ben John Fisher," Members, ASCE, 
and Stuart M., ASCE 


INTRODUCTION 


One the dominant features many tidal inlets the ebb tidal shoal. 
Generally characterized horseshoe shape, this lobe sediment often controls 
and turn controlled the tidal currents and channels (5). The ebb tidal 
shoal the zone deposition for sediment carried tidal currents from 
the interior bay estuary and longshore currents. The material moved 
suspension and bedload the ebb discharge deposited the flow 
slows with increased depth, lateral spread, and mixing. general, the size, 
geometry, and position the shoal appear depend the texture its 
sediment, strength tidal and longshore currents, inlet morphology, and 
nearshore wave climate. 

While significant work has been done the morphology ebb tidal shoals 
(5,2,6), little has been published the relationship this deposit the controlling 
hydrodynamics. Ozsoy (7), Dean and Walton (4), and others have suggested 
both theoretical and experimental approaches, and those studies guided our 
own efforts. Many ebb tidal flows can characterized turbulent jets (9,7) 
which can analyzed relatively straightforward manner, both theoretically 
and the laboratory. investigated the development the ebb tidal shoal 
the sediments which were deposited after transport such jet. 

The mechanics shoal formation profound interest coastal engineers. 
The manipulation inlet geometries and dredging navigation channels are 
but two many design problems which either directly indirectly involve 
the ebb tidal shoal. order address these questions adequately, predictive 
model for shoal formation needed. Our study intended basic step 
the development this model. Using small laboratory basin, ebb shoals 
were formed with steady jet. Waves, reversing flows (flood tide), and sloping 
bottom were not studied, our conclusions cannot applied directly actual 
inlets. Rather, this initial study (8) was designed evaluate the appropriateness 

‘Assoc. Prof., Dept. Civ. Engrg., Clemson Univ., Clemson, S.C. 

Asst., Dept. Civ. Engrg., Clemson Univ., S.C. 

Note.—Discussion open until April 1982. extenu che closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication August 22, 1980. 
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the jet hypothesis the laboratory model. Positive results from these initial 
experiments will provide foundation for more complete laboratory study. 


Test 


The 3.7 3.7 0.6 model basin used this investigation shown 
schematic diagram Fig. Water was withdrawn through perforated 
weir that served minimize undesirable recirculation currents. The pump 
discharge flowed around the basin through valve and flowmeter and then 
through gravel filters which eliminated ripples and waves and allowed smooth 
flow through the inlet. The width the inlet was adjusted changing the 
separation two 1.9 thick partitions. Test inlet widths varied between 
1.3 and 5.1 cm, water depths varied from 4.3 15.2 cm, and flow 

After the inlet geometry was adjusted, sediment was placed directly front 
the inlet one two ways. one method, mound (175 cm’) 
front the inlet. the other method, the sediment was placed 0.3 thick 
over 0.6 wide 0.9 long area front the inlet. Duplicate tests 
showed that neither the amount sediment nor method placement affected 
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FIG. 1.—Schematic Diagram Model Basin Hydraulics Showing Developing Ebb-Tidal 
Shoal 
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the results. During test, the valve was adjusted the desired flow rate, 
and the shoal was allowed grow until length change greater than 


TABLE 1.—Summary Shoal Sediment Parameters 


Critical 
velocity, 
milli- meters 
Sediment type Specific meters 
(test numbers) Gravity 
(1) (5) 
Sand (1-18) 
Sand (19) 
Crushed walnut shells 
G-2 crushed glass (27-28) 
Crushed brick (29-30) 


FIG. 2.—Locations Important Dimensions Shoals 


one percent one hour was achieved. The shape and size the shoal (Fig. 
was then photographed diagrammed. 
During the course the test program, total thirty shoals were formed 


tir 
— 
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from five sediment types various particle sizes and specific gravities (Table 
1). Steady-state shoal dimensions for the thirty tests are given Table 


TABLE 2.—Test Results 


Flow 
rate, 
cubic 

centi- Outside Outside Inside Inside 
Inlet Water meters length, width, length, width, 

centi- milli- centi- centi- centi- centi- 
meters meters meters meters meters meters 


(2) (3) (5) (6) (7) (8) 


21.08 
6.86 
8.13 

15.75 

22.61 
6.10 

14.73 

18.29 

15.24 

13.46 

27.69 
2.79 
9.40 

14.22 

16.00 

23.37 

24.13 
7.11 

19.81 
6.35 
7.87 

20.32 
8.38 

24.64 

14.22 

2.29 
12.45 
19.30 
17.78 
17.78 


Regardless the type sediment test situation, the classical horseshoe 
shape prevailed throughout the shoal formation. 


For the thirty tests conducted, shoal dimensions varied substantially 
sediment type, discharge velocity, u,, inlet width, W,, and depth, 


) 
2.39 10.16 1.04 76.20 26.67 69.85 
2.54 10.16 29.46 17.02 16.76 
2.54 10.16 30.99 15.24 22.00 
2.54 10.16 59.44 23.62 42.42 
2.54 10.16 1.20 85.09 28.45 77.72 
2.54 12.70 23.37 11.68 17.27 
2.54 12.70 1.06 45.21 22.10 37.85 
2.54 12.70 1.18 57.40 25.40 51.05 
3.17 12.70 1.14 47.24 22.10 38.61 
3.17 12.70 1.14 50.55 20.32 41.15 
3.17 8.89 1.20 99.57 33.78 92.56 
3.81 8.89 11.43 7.62 3.56 
3.81 8.89 33.78 16.51 20.83 
3.81 8.89 1.20 47.24 22.86 34.29 
3.81 8.89 1.20 45.47 23.37 39.88 
3.81 8.89 1.20 85.09 29.46 81.28 
3.81 8.89 1.20 85.09 32.51 74.68 
15.24 1.20 22.61 14.73 15.49 
3.81 12.70 1.23 53.85 28.45 47.24 
1.65 6.55 .26 43.43 11.94 32.26 
2.46 35.31 14.99 26.92 
2.59 10.36 73.66 27.43 60.20 
1.73 6.93 46.48 14.22 39.37 
3.45 13.84 1.06 77.47 30.48 72.14 
2.31 9.25 76.20 21.08 68.58 
1.29 4.37 18.54 6.86 11.18 
2.54 12.37 45.97 19.05 40.39 
3.81 12.70 1.20 50.29 26.42 44.70 
3.17 12.04 1.07 46.74 24.38 40.13 
2.54 11.73 1.14 67.31 21.59 59.18 
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500 sand 
walnut 
glass 
crushed brick 


FIG. 3.—Nondimensional Maximum Shoal Width Versus Nondimensional Initial 
Average Velocity for all Sediment Types 


500 sand 
420 sand 
glass 

crushed brick 


FIG. 4.—Nondimensional Maximum Shoal Length Versus Nondimensional Initial 
Average Velocity for all Sediment Types 
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performing dimensional analysis the pertinent parameters, was necessary 
relate the sediment type one more geometric hydrodynamic parameters. 

Based classical bedload transport studies, the Shields entrainment function, 
relatively constant for the range particle Reynolds numbers, u,d/v, 
tested. Thus 


2 
u, 


constant 


shear velocity; and the particle specific gravity. For fixed bed roughness, 
proportional the average fluid velocity mid-depth, (measured 
electromagnetic current meter for this study), thus 


u 


50 


which corresponds the fluid velocity sufficient initiate particle motion. 
simple procedure determine whether the shoals formed the present 
study satisfy Eq. The average fluid velocity above the terminal lobe 
shoal equilibrium was measured and used with the sediment properties 
calculate the constant Eq. (Table 1). While there some variation, 
not large and the implication that for the initial analysis sediment properties 
(of diameter and specific gravity) can replaced single variable, 


W,, 


Employing the Buckingham Theorem, dimensional analysis Eq. 


gives, after observing from test results that water depth had little effect 
shoal size 


which either shoal length width. 

Nondimensional plots shoal dimensions are given Figs. and The 
ability successfully normalize results from all tests opens the possibility 
that studies such these can used assist formulating predictive models. 

the extent that the laboratory shoals are similar shape those 
actual inlets, several observations can made with regard their response 
the controlling parameters. The slope both Figs. and sufficiently 
small indicate relatively close correlation shoal size and inlet velocity. 
For example, doubling the velocity ratio from 1.2 2.4 triples shoal 
length (from 30). Moreover, for the range test conditions, 
the equilibrium shoal dimensions vary linearly with inlet velocity. 


Having developed actual laboratory shoals using conditions usually satisfied 
one-dimensional jet analysis, analytical and experimental results were 
compared. Since sediment movement related the bottom shear stress, which 
turn proportional the jet velocity, was assumed that isopleths 
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shear stress which would likely correspond shoal position could approximat- 
jet isotachs. These isotachs were obtained from classical 1-D jet theory, 
based the governing equation control volume form. 

Conservation streamwise (x) momentum written 


which b(x) the jet half width; the velocity the direction; 
and the lateral coordinate. customary assume that velocity profiles 


FIG. 5.—Critical Velocity Shown with Shoal for Test Plotted 0.3 meter 
Square Grid 


exhibit dynamic similarity—our data also indicated this—or 


which u,, the maximum centerline velocity. Substitution Eq. into 


which the constant 


(5) 
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Our tests agreed with previous results (1) that increases linearly with 


x as 

which 0.22; and velocity profiles could approximated cubic 


w 
N 


FIG. 6.—Test Results Showing that Velocity above Shoal (at mid depth) not 
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overall check the conformity our inlet tests with accepted 
jet theory, Eq. was substituted Eq. u,, replaced the critical velocity 
the end the shoal u., X,,, the shoal length, and was calculated 
from 


The average value computed this way from the first tests was 0.32, 
opposed theoretical value 0.37 for Eq. 10. summary, the jet discharges 
produced during this study exhibited growth rates, dynamic similarity velocity 
profiles, and values reasonable agreement with that generally accepted 
for 1-D jets. 

After demonstrating that the inlet jets produced this study were basically 
one-dimensional, calculation isotach shapes was straightforward. Eq. was 
used with and combined with and obtain isotach locations. 
Fig. shows results for typical sand shoal. Results this exercise indicated 
that all the laboratory shoals were much wider than the isotach corresponding 
the critical velocity taken above the terminal lobe. fact, the theoretical 
length-to-width ratios varied from 7.4 14.6, while the experimental ratios 
varied from 1.3 5.1 and averaged one-fifth the theoretical values. Velocity 
measurements above one the shoals (Fig. substantiate this result—the 


shoal obviously was not isotach. Ozsoy (7), using two-dimensional jet theory 
combined with sediment transport model, predicted length width ratios 
about 10, also much greater than those observed this study. The marked 
difference indicates possibility different transport mechanisms for the terminal 
lobe and the channel margin linear bars, result also suggested Butakov 


(3). 


All shoals generated during this study contained the characteristics generally 
attributed those occurring tide-dominated inlets. This was true regardless 
the flow rate, inlet width, sediment type used. The dimensions the 
laboratory shoals were functions inlet width, initial average velocity, and 
sediment type. 

For this initial test program, the inlet width and sediment 
were successfully used normalize the relation between shoal size and discharge 
velocity. The resulting plots imply that some the dominant parameters have 
been identified. The success analysis encouraging and suggests that 
further research may permit the development truly predictive model. 

obtain better understanding the growth process, several time-lapse 
movies were made shoals from discharge initiation their equilibrium position. 
Observations showed that the length width ratio (aspect ratio) shoals 
their early growth phase was larger than the equilibrium value. For example, 
the aspect ratio one shoal varied from 5.7 the early stages its development 
3.7 its final equilbrium position. While these early aspect ratios were 
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higher than the previously mentioned equilibrium values, they were still 
not high the theoretical values obtained from the isotach shape. More 
detailed studies are needed determine there fact substantial difference 
transport mechanisms the channel margin linear bars and the terminal 
lobe. 


Results the test program and analysis can summarized follows: 


Equilibrium shoal length and width increase direct proportion the 
average inlet velocity. 

Inlet width and sediment velocity measured mid-depth are 
appropriate scale factors for the laboratory shoals. 

Two-dimensional theoretical isotach shape does not satisfactorily predict 
the equilibrium shoal shape. 


More definitive conclusions await the generation larger laboratory shoals 
which use movable bed and include longshore currents and sloping bottom. 
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Notes 


provide place witain ASCE for publication technical ideas that have not advanced, 
yet, the point where they warrant publication Proceedings paper Journal, the 
publication Technical Notes was authorized the Board Direction October 16-18, 
under the following guidelines: 


original manuscript and two copies are submitted the Manager Technical 
and Professional Publications, ASCE, 345 East 47th Street, New York, N.Y., 10017, along 
with request the author that considered Technical Note. 

The two copies will sent appropriate Technical Division Council for review. 

the Division Council approves the contribution for publication, shall returned 
Society Headquarters with comments. 

The technical publications staff will prepare the material for use the earliest possible 
issue the Journal, after proper coordination with the author. 

Each Technical Note not exceed pages the Journal. approximation, 
each full manuscript page text, tables, figures the equivalent one-half Journal 
page. 

The Technical Notes will grouped special section each Journal. 

Information retrieval abstracts and key words will unnecessary for Technical Notes. 

The final date which Discussion should reach the Society given footnote 
with each Technical Note. 

Technical Notes will not included 7ransactions. 

10. Technical Notes will included ASCE’s annual and cumulative subject and author 
indexes. 


The manuscripts for Technical Notes must meet the following requirements: 


Titles must have length not exceeding characters and spaces. 

The author’s full name, Society membership grade, and footnote reference stating present 
employment must appear the first page the manuscript. Authors need not Society 
members. 

The manuscript submitted original copy (with two duplicates) that typed 
double-spaced one side 8-1 /2-in. (220-mm) (280-mm) white bond paper. 

All mathematics must typewritten and special symbols must properly identified. 
The letter symbols used must defined where they first appear, figures text, and arranged 
alphabetically Appendix.—Notation. 

Standard definitions and symbols must used. Reference must made the lists 
published the American National Standards Institute and the Authors’ Guide the 
Publications ASCE. 

Tables must typed double-spaced (an original ribbon copy and two duplicate copies) 
one side 8-1/2-in. (220-mm) (280-mm) paper. explanation each table 
must appear the text. 

Figures must drawn black ink one side (220-mm) (280-mm) 
paper. Because figures will reproduced with width between in. (76 mm) 4-1/2 
in. (110 mm), the lettering must large enough legible this width. Photographs must 
submitted glossy prints. Explanations and descriptions must made within the text 
for each figure. 

References cited text must typed the end the Technical Note alphabetical 
order Appendix.—References. 

Dual units, i.e., U.S. Customary followed (International System) units parentheses, 
should used throughout the paper. 


TECHNICAL NOTES 


SMOOTHING PREVENT 
NUMERICAL INSTABILITY 


INTRODUCTION 


Numerical simulation sediment transport and the changes plan form 
the shoreline near coastal structures may produce results where the bottom 
slope exceeds the maximum slope which would expected occur the 
real (prototype) shoreline. This may cause instability the numerical solution 
and prevent accurate simulation the transport processes. numerical 
subroutine has been developed redistribute the sediment points where 
the numerically produced bottom slope exceeds the maximum expected slope. 
This subroutine has dependence time. The smoothing accomplished 
numerically after the completion one more iterations time-dependent, 
numerical simulation physical processes, e.g., sediment transport associated 
with wave action currents. should noted that the smoothing subroutine 
will smooth any input data reach maximum slope between points, even 
random numbers. This subroutine should not used substitute for good 
simulation the physical processes, but may used smooth errors resulting 
from the inexact nature numerical modeling. 


assumed here that rectangular array used for the numerical computation. 
The results presented can modified for cases where rectangular array 
not used, e.g., near boundary. nine-point rectangular array used for 
the computations each grid point (Fig. 1). The water depths are expressed 
positive numbers for ease computation. 

The first step the computation the determination potential sediment 
movement, between point (x,, y,) the center the array and each 
the surrounding points. difference depth, which would result maximum 
slope, S,,, expected occur between the two points for the known physical 
conditions, compared with the difference depth produced the numerical 
simulation the physical processes moving the sediment. This gives 

Engr., U.S. Army Coastal Engrg. Research Center, Kingman Bldg., Ft. Belvoir, 
Va. 22060. 
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r, = Ay, at 


r,= Ax, at (X41> Y,) 


cal Computation 


which the water depth the point 
another point the array, and r,, the distance between the two points, given 


Where sediment motion assumed occur, and set equal 


Water Surfoce 


FIG. 1.—Nine-Point Rectangular Array and Representation Sea Bottom for Numeri- 


zero those points. 
where motion would occur, the direction potential motion 
and from the remaining points must determined. This shown 


which converts negative number, v,, where D,, i.e., where the 
potential sediment movement away from the center the array. 

then necessary determine the interrelationships the potential sediment 
motions. This done first expressing the potential sediment motions 
terms their relative influence 
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(2) 
(3) 


After considering all surrounding points and setting 


TECHNICAL NOTES 


Next, the total sums the positive and negative potential sediment movements, 
and are expressed terms equivalent numbers maximum poten- 


Equilibrium 
Position 


FIG. 2.—Schematic Depth Correction, 


v 


Equilibrium 


FIG. 3.—Schematic Depth Correction 
tials, and Tespectively, that 


and 
(v; ) max 
determine correction coefficients, and c,, correct the positive and 
negative potential sediment movements, and respectively, the case 


(increase center) (decrease exterior points) 
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unity. The equation for then developed noting increase (on left) in- 
crease (at center) decrease (on right) 


2 max 
. . (8; 


Noting from Fig. that 


and substituting the right-hand side from equation 


n 
2 8) max 
| max 
P 


Using the appropriate values and C,, now possible correct 
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the corrected depth. For cases which D,> and8,>0, from Eq. 


(remembering that has negative value). From conservation mass 


Coerricients 


Because the numerical computation inexact, necessary consider 
how the corrections Eqs. 20-22 affect the data. Using the full correction 
shown the equations might, some instances, cause 
the input data; i.e., too much sediment moved the numerical subroutine 
given point because the point-by-point computation through the matrix. 
This may cause excessive flattening the bottom slope (to slope less 
than which will not corrected successive iterations. This could 


partly compensated for using smaller correction coefficients, and 
given 


which 1.0, and using larger number iterations through the matrix. 
This was tested using matrix. Use the full-correction coefficients 
required iterations through the matrix reach solution. Where 0.6, 
iterations were required, and where 0.3, iterations were required. 
While use the smaller correction coefficients provided some improvement 
the final values water depth individual points, this was offset 
substantial percentage increase computer time. The actual cost increase 
computer time depends the complexity the problem where the subroutine 
used, and the number times the subroutine called. The desirability 
using smaller coefficients must determined the individual programer. 
smaller coefficients are used, the necessary condition 


has been assumed that sediment motion will occur points which 
However, attempt made here account for this additional sediment motion. 
Direct computation each grid point the network would require significant 
increase computational time, i.e., recomputing depth differences and readjust- 
ing depths through number iterations, and the smoothing computations 
from point-to-point plus successive iterations through the entire network should 
account for the resulting motion. 
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DISCUSSION 


Closure Subrata Chakrabarti,* ASCE 


The discussions Ramberg and Christodoulou are appreciated. 

There seems some confusion regarding the data reduction procedure 
adopted the paper. test record typically was ten cycles long. Irregular 
stream function theory was used predict the water particle kinematics from 
the wave profile. Each cycle the profile was fitted with the theory and 
the water particle velocity and acceleration profiles were derived from it. Since 
there was some variation the wave profile from one cycle the next, the 
values the kinematics, e.g., u,, were different from cycle cycle. However, 
the variation these values was generally small and account for the clustering 
the data points Fig. These variations the C,, and values within 
one wave record were much greater than any errors the measurement. 
the absence the measured kinematics the flow field near the cylinder 
difficult attribute the cause the variations physical factor. However, 
possible that the flow field the vicinity the cylinder constantly 
changing due change the vortex shedding process from the cycle cycle 
variation the free surface profile. This, turn, alters the values C,, and 
controlled planar harmonic flow (e.g. Sarpkaya’s experiment) such 
cycle cycle variation avoided. 

Regarding comparison the hydrodynamic coefficients from the wave tank 
test with those from the planar flow results, should noted that only the 
mean values data were used. Referring Fig. may seen that the 
spread the data, particularly the values, large. the lower end 
the envelope, the values are appreciably lower than Sarpkaya’s results, 
while the upper end values are higher. Ramberg’s (13) comparison was made 
only with eight data points and cannot called conclusive. Moreover, his 
data were averaged over the length the vertical cylinder (since total forces 
were measured) and his discussion applied range values where 
the drag force contribution quite small compared with the inertia force 
contribution and, such, amenable substantial error. For example, 
the drag force about 10% the total force, error the measurement 
may produce large error the computed drag coefficient. The writer agrees 
with Ramberg’s explanation the variation the two results partly 
contributing the three-dimensionality the flow field and different kinematic 
gradients due orbital variation. 

should further noted that experiment with vertical cylinder, 
the coefficients should the measured force data from small segment 
the cylinder. Any conclusion derived from averaged values over the length 
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the cylinder influenced the averaging the complex three-dimensional 
flow field. 

Regarding the discussion Christodoulou, is, indeed, difficult reach 
any fruitful conclusions regarding the influence C,, and from the 
present study. This primarily due the small range While the correlation 
between the results from the wave-flow here and the two-dimensional harmonic 
flow for the appropriate values may termed fair, the difference 
still large, especially view the scatter the C,, and data. more 
controlled experiment wave tank with larger variation needed 
for the verification these results. This, however, extremely difficult 
because the limitation wave generation the test facility. the 
other hand, field test data are difficult reliably obtain and their results difficult 
interpret, due additional environmental effects the test cylinder. 

The writer agrees with Christodoulou that the correlation the total forces 
using data from the local forces the same test cylinders somewhat biased. 
Nonetheless, establishes the accuracy the measuring instruments since 
two types instruments were used measuring the total and local forces 
and the local forces are much more delicate measure. Since more than one 
test cylinder and test set-up was used, provides greater confidence the 
data. also determines the appropriateness and usefulness the mean curves 
for C,, and through large scatter data. The correlation forces for 
the 1.5 in. diameter tube (Fig. poorer, not because contributed fewer 
points, but because the computation applied the mean curves the higher 
range where data are sparse and subject larger error. would 


worthwhile correlate the data using results derived from other independent 
tests. 


Errata 


The following corrections should made the original paper: 


“November, 1980, Liu Teh-fu and Feng-shi (Proc. Paper 15862). 
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